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We observe experimentally that ensembles of quantum dots in three-dimensional (3D) photonic 
crystals reveal strongly nonexponential time-resolved emission. These complex emission decay curves 
are analyzed with a continuous distribution of decay rates. The log-normal distribution describes 
the decays well for all studied lattice parameters. The distribution width is identified with variations 
of the radiative emission rates of quantum dots with various positions and dipole orientations in the 
unit cell. We find a striking sixfold change of the width of the distribution by varying the lattice 
parameter. This interpretation qualitatively agrees with the calculations of the 3D projected local 
density of states. We therefore conclude that fiuorescence decay of ensembles of quantum dots is 
highly nonexponential to an extent that is controlled by photonic crystals. 



I. INTRODUCTION 

Control over spontaneous emission from ensembles of 
excited light sources is of great interest for many applica- 
tions, such as miniature lasers, light-emitting diodes, ^'^ 
and solar cellsi^ The rate of spontaneous emission is de- 
termined not only by the internal nature of emitters 
but also by their environment 4*^ According to Fermi's 
golden rule, this rate is proportional to the local ra- 
diative density of optical states (LDOS), to which the 
emitters couplej^ii This projected LDOS counts, at given 
frequency and orientation of the transition dipoles, the 
number of electromagnetic states at the locations of the 
emitters. It has been predicted that periodic dielectric 
structures, known as photonic crystals, can be used to 
radically change the LDOS.^ The main research goal is 
the achievement of a photonic band gap, i.e., a range 
of frequencies where no electromagnetic states exist in- 
side the crystal, irrespective of the location. It has also 
been predicted that a much weaker requirement than a 
gap suffices to suppress spontaneous emission>& by plac- 
ing sources at judicious locations in the crystal unit 
cell where the LDOS vanishes. Since the frequency- 
integrated number of states is conserved, one expects the 
LDOS to be strongly increased at some frequencies out- 
side such a pseudogapi^'^'^ d°di This means that photonic 
crystals may completely control the emission rates be- 
tween complete inhibition and strong enhancement even 
in the absence of a photonic bandgap. Since photonic 
crystals are a natural platform for solid-state emitters 
such as quantum dots^^ii^ii^ii^ such control of sponta- 
neous emission is relevant to applications in quantum in- 
formation. 

Most theoretical papers on spontaneous emission in 
photonic crystals concern single light sources i^i^ d^i^^ In 
the case of a weak emitter-field interaction, one expects 
to see a single-exponential decay curve with a slope equal 
to the decay rate or inverse lifetime. However, experi- 



ments on ensembles of emitters often show nonexponen- 
tial decays. Such complex decay dynamics can be due 
to four reasons: (i) Emitters experience different LDOSs 
when they are distributed over different positions and 
dipole orientations in the unit cell of a photonic crys- 
tal, (ii) It has been predicted that single sources reveal 
nonexponential decay due to van Hove singularities in 
the LD0S,-2 (iii) Nonexponential decay may appear if 
the emitters have more internal levels than the usually 
considered two levels, (iv) Temporal fluctuations of the 
emitters' environment on time-scales larger than the flu- 
orescence lifetime can lead to apparent nonexponential 
decays. At any rate, it is an open challenge to interpret 
complex nonexponential decay curves, in particular, to 
obtain information on the local density of states. 

In this study, we investigate time-resolved spontaneous 
emission from an ensemble of light sources distributed 
over a well-defined set of positions in the unit cell. We in- 
terpret the emission data with a continuous distribution 
of emission rates. This distribution is identified with the 
distribution of the LDOS over all positions r; sampled by 
sources with fixed emission frequency and random dipole 
orientations. The distribution width shows a striking six- 
fold variation with the varying crystal lattice parameter, 
in qualitative agreement with intricate calculations of the 
three-dimensional (3D) LDOS. This study opens an av- 
enue to the analysis of time-resolved emission from en- 
sembles of light sources in complex photonic systems. 

II. EXPERIMENTAL DETAILS 

We have studied room-temperature spontaneous emis- 
sion from sources embedded in fee inverse opals consist- 
ing of air spheres in a titania (Ti02) backbone shown 
in Fig. [TJa) . These 3D photonic crystals strongly inter- 
act with light . ^^d"^ Extensive details on the fabrication 
and characterization of the inverse opals are reported in 
Ref. [H. The lattice parameter a is determined from the 
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FIG. 1: (a) Electron-microscope image of the (111) surface 
of an inverse opal consisting of air spheres in Ti02 with lat- 
tice parameter a = 425 nm. (b) Schematic on quantum dots 
(black spots) on the internal titania-air surfaces at symmetry- 
inequivalent positions r; in the crystal unit cell. 



measured central wavelength of the lowest stop band Ac: 
a — V3Ac/(2n), where n = 1.27 ± 0.15 is the average 
refractive index consistent with 10 to 20 vol % Ti02. We 
have studied 15 different samples with lattice parameters 
ranging from a — 255 ± 10 nm to a = 760 ± 20 nm. 

We use CdSe-ZnSe (core-shell) colloidal quantum dots 
(QDs) as light sources because of their high fluorescence 
quantum efficiency and narrow homogeneous spectral 
width; The emission spectrum of the QDs is centered 
at A = 610 nm, which is determined by the average di- 
ameter of the nanocrystals of 4.5 nm. The process of the 
liquid infiltration of the photonic crystals with the QDs 
is described in Refs [l2fl7l : the QDs precipitate at posi- 
tions ri that are random but within a well defined set, on 
the internal surfaces of the air spheres inside the inverse 
opals (see Fig.[T]D), with an estimated low density of four 
QDs per air sphere. 

The QDs are excited at A = 447 nm with a diode laser 
(Picoquant) emitting 90 ps pulses with 20 nJ/pulse^ii 
We record fiuorescence decay curves of QD emission 
with a microchannel plate photomultiplier tube detector 
(Hamamatsu R3809U) using the time-correlated single- 
photon counting method. The decay curves are his- 
tograms of the arrival time of a photon emitted after 
the laser pulse, obtained over many excitation-emission 
cycles with a resolution better than 100 ps. The slope 
of the decay curves yields a decay rate 7 = jrad + Inrad, 
which is the sum of the spontaneous-emission decay rate 
7rad and the nonradiative rate "fnrad depopulating the 
excited states of the QDs. 



III. RESULTS AND DISCUSSION 

In Fig. Owe show time-resolved spontaneous emission 
from the QDs in inverse opals with three different lat- 
tice parameters a, i.e., different reduced frequencies a/ A. 
The data were collected at A = 615 nm within a narrow 
range AA = 3 nm to select the same population of QDs 
with identical emission properties on each sample. For 
inverse opals with a = 255 nm, A = 615 nm is in the 
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FIG. 2: Fluorescence decay curves recorded at A = 615 nm 
and T = 295 K from QDs in inverse opals with lattice pa- 
rameters of a = 425 nm (1), a = 255 nm (2), and a = 540 
nm (3); and from QDs in a different chemical environment, a 
chloroform suspension (S). The solid lines are fits of the log- 
normal distribution of decay rates to the data. The goodness 
of fit Xred varies from 1.1 to 1.4, close to the ideal value of 1. 
The dashed line in (3) shows a heuristic biexponential model, 
which does not agree with the data (Xred > 1-9)- 

low- frequency limit, where the frequency dependence of 
the LDOS is known to show an a;^-behavior.S' — We see 
that the spontaneous emission in a sample with a — 425 
nm is inhibited compared to the reference. Conversely, in 
a sample with a — 540 nm, the decay rate is enhanced. 
As in any photonic-crystal environment, the backbone 
also fluoresces, which here distorts the signal at short 
times. We have carefully removed the Ti02 signal from 
the measured decay curves since we know its spectrum 
and decay curve from the measurements on an undoped 
inverse opal; the backbone has a count rate less than 12 % 
of the QD signal. We exclude the possibility that QDs at 
the sample surface contribute any measurable signal: the 
analysis of separate angle-resolved measurements reveals 
an excellent agreement with theory for sources emitting 
in the bulk of photonic crystalsfi^ QDs on the surface are 
rinsed off after infiltration. Since we also verified that the 
pump beam was not Bragg diffracted, the pump intensity 
is maximum inside the samples at a distance of several 
Bragg attenuation lengths from the surface due to light 
diffusion. Therefore, Fig. [2] demonstrates time-resolved 
emission from an ensemble of QDs controlled by photonic 
crystals, over a much larger time span than in previous 
experiments J^ii^ii^ii^i^ 

A remarkable feature in Fig. [2] is that the decay curves 
from the QDs in the inverse opals strongly deviate from 
a single-exponential decay. To explain this observation, 
we consider the four reasons discussed above: i) Since 
the QDs are distributed over positions and dipole orien- 
tations in the unit cell (see Fig. [1)3) , they should indeed 
experience different LDOSs. ii) Observation of nonexpo- 
nential decay due to van Hove singularities in the LDOS 
requires single-dot experiments, which is not the case 
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here, iii) Even though the colloidal QDs are not true two- 
level systems, their emission decay is close to being single 
exponential,— as confirmed in Fig. [21 iv) It was suggested 
in Ref. ^23 that temporal fluctuations of the environment 
surrounding the QDs induce a distribution of nonradia- 
tive decay channels. In our experiments, however, the 
nonradiative rates ^nrad hardly vary from sample to sam- 
ple because QDs from the same batch are used, and the 
photonic crystals are chemically identical. We observe 
only minute differences of decay-curve slopes among the 
samples with the same lattice parameters, which indi- 
cates that the temporal fluctuations are identical for all 
samples. Therefore we can safely attribute the observed 
variations of the nonexponential decay curves to a dis- 
tribution of radiative decay rates "/rad as a result of a 
spatial and orientational variation of the LDOS. 

To interpret the complex, nonexponential decay 
curves, we propose a different line of attack by model- 
ing the curves with a continuous distribution of decay 
rates: 

poo 

lit) - /(O) / <^(7)e-^*d7, (1) 

where (/)(7) is a distribution of decay rates with dimension 
of time. The fluorescence intensity /(t) is proportional to 
the time-derivative of the concentration of excited emit- 
ters. Therefore, describes a distribution of the con- 
centration of emitters with a certain 7, weighted by the 
corresponding jrad-— This approach has two advantages: 
first, it enables us to explain intrinsically nonexponen- 
tial decay curves, and second, the distribution containing 
physical information on decay rates is readily available, 
which is essential when treating an ensemble of emitters. 
We use the log-normal distribution function 

0(7)=Acxp(-i^^!i2f^), (2) 

where 7mf is the most-frequent decay rate correspond- 
ing to the maximum of 0(7), u; is a dimensionless width 
parameter that determines the distribution width at 1/e: 

A7 = 27j\/F sinhiu. (3) 

A is the normalization constant, so that 0(7)^7 — 1. 
The important features of the log-normal distribution 
are that the logarithmic form of the distribution func- 
tion excludes unphysical negative decay rates and that 
it is specified in terms of only two free parameters, jmf 
and A7. Other multiexponential models are the heuristic 
biexponential model and the Kohlrausch stretched ex- 
ponential model that has been employed to QDs out- 
side photonic crystals,^^^^^ Figure O shows that the bi- 
exponential model does not match our data, even though 
more free parameters are involved. The stretched expo- 
nential model does not match our data eitherf2^ which 
again confirms that the variations we observe are due to 
LDOS effects in photonic crystals and not to complex 




Decay rate Y = Y^^^+ Ynrad ["s"'] 

FIG. 3: Decay-rate distributions </>(7) for the inverse opals 
with lattice parameters a — 425 nm (1), a = 255 nm (2), and 
a — 540 nm (3), corresponding to the data shown in Fig. [2] 
The inset shows 0(7) for a — 255 nm in the semilogarithmic 
scale. Clear modifications of A7 and 7a/f with varying lattice 
parameter of the inverse opals are seen. 

emission properties of the QDs. In Fig. [5] it is seen that 
the log-normal distribution model (solid curves) provides 
an excellent description of the experimental data. 

Figure [3] shows the resulting decay-rate distributions 
for three lattice parameters. It is remarkable that the 
log-normal distribution model provides an excellent ex- 
planation for all reduced frequencies a/X studied, which 
will be seen below to agree with calculations. Compared 
to the low- frequency reference (a = 255 nm), the max- 
imum of the distribution (/'(7) is shifted to lower decay 
rates for the crystal with a = 425 nm and to higher rates 
for the crystal with a = 540 nm. These shifts are a clear 
demonstration of a photonic effect of the inverse opals 
on the ensemble of embedded emitters. In Fig. [31 we see 
a dramatic change of the width A7 of the distribution. 
The large width of each distribution is identified with 
the variation of the radiative emission rates due to ori- 
entational and positional ri dependences of the LDOS at 
each lattice parameter. Consequently, the decay rates of 
individual QDs are much more strongly modified by the 
photonic crystal than the most-frequent rate 7a/f of the 
ensemble. 

In Figs [3^ and [Hd we have plotted the resulting values 
for A7 and jmp versus the reduced frequency. Let us 
briefly consider jmf- because the nonradiative part of 
the total decay rates does not change with the lattice pa- 
rameter, the change of 7mf - Imp purely radiative and 
related to an averaged photonic-crystal LDOS. We com- 
pare the experimental data to the calculated density of 
states (DOS) - the unit-cell average of the LDOS (dashed 
curve). The measured variation of ^mf is seemingly in 
good agreement with the DOS. Both inhibited and en- 
hanced decay rates are observed, and the experimental 
variation in 7mf amounts to a factor of 3. In contrast, 
the width A7 shows a striking sixfold variation, much 
larger than the change of 7AfF (Fig- W^)- Already in the 
low-frequency limit, a/A = 0.4, there is a spatial varia- 
tion of the radiative rate (A7 > 0) because the QDs dis- 
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FIG. 4: (a) Width of the decay-rate distribution A7 vs. 
lattice parameter at a fixed emission wavelength (triangles, 
dashed curve is a guide to the eye), (b) Difference between 
7A/F measured on a photonic sample and 'y^iF ~ 0-05 ns"'^ of 
the low- frequency reference (squares). The error bars are esti- 
mated from the largest difference between data on the samples 
with similar a. The dashed curve represents the relative DOS 
in the inverse opals. The relative (L)DOS in the inverse opals 
is the (L)DOS divided by DOS in a homogeneous medium 
with the same average refractive index (n = 1.27). (c) Rela- 
tive LDOS at two positions on the internal Ti02-air surfaces 
projected on dipole orientations parallel (||) or normal (±) to 
the internal surface; at ra = 0.144(1,1,2) and rt = 0.2(1,1,1). 



725 reciprocal-lattice vectors in the iJ-field plane-wave 
expansion method i^"'^^ The inverse opals were mod- 
eled as close-packed air spheres surrounded by shells of 
Ti02 (e = 6.5) with cylindrical windows between neigh- 
boring spheres. The model agrees with prior optical 
experimentsi^ The integration over wave vectors k was 
performed by representing the full Brillouin zone by an 
equidistant k-point grid consisting of 291 416 points. 
The results in Fig. [4j: reveal a strong dependence of the 
LDOS both on the position in the crystal unit cell (com- 
pare curves 1 and 2) and on the dipole orientation (com- 
pare curves 1 and 3). It is remarkable that in the rele- 
vant a/A range, the dependence of the LDOS on the lat- 
tice parameter is the same at both positions r; and both 
orientations, and even as the unit-cell averaged LDOS 
(Fig. |3)d). This result agrees with the observation that 
all measured decay curves are successfully modeled with 
the same log-normal shape of the decay-rate distribution. 
Because the LDOS for dipoles perpendicular to the inter- 
face is inhibited and nearly constant at all reduced fre- 
quencies, whereas the LDOS for parallel dipoles strongly 
varies, we propose that the width of the LDOS distri- 
bution has a similar frequency dependence as the LDOS 
itself. This notion agrees with our observation that A7 
tracks the behavior of 7m_f- A quantitative comparison 
of our data to the calculated LDOS is a challenge, since 
detailed knowledge is needed on the relation between ^rad 
and 7 to infer the true radiative decay-rate distribution 
(see Ref. Hsl). Qualitatively, the calculated LDOS reflects 
the main features of our experiments. 



IV. SUMMARY 



tributed over inequivalent positions in the unit cell cou- 
ple to different electric fieldsi^ At the frequencies of the 
i-gap, a/A « 0.7, the radiative rate ^rad is inhibited in 
most places in the unit cell that are occupied by the QDs, 
as confirmed by a low continuous-wave (cw) count rate 
of only w 2.5 kHz. Therefore, the observed narrow width 
A7 = 0.1 ns~^ is a measure of the distribution width 
of nonradiative rates. In contrast, at a/A = 0.88, A7 is 
strongly increased; here the cw count rate is w 56 kHz at 
similar experimental conditions, in agreement with an en- 
hanced ^rad- We therefore conclude that the large widths 
A7 are determined by a broad distribution of radiative 
emission rates IS.'yrad that are proportional to a broad 
distribution of the projected LDOS at fixed frequency. 
Hence, the width A7 is a much more characteristic pa- 
rameter to describe the ensemble emission in 3D photonic 
crystals than the usually used average ratesjiS, 

We have managed to perform intensive computations 
of the 3D LDOS at two representative positions in 
the unit cell at the Ti02-air interface (see Fig. Wp)- 
The LDOS shown in Fig. Uj: was calculated for dipole 
orientations parallel and perpendicular to the Ti02- 
air interface. The calculations were performed using 



We have successfully explained highly nonexponential 
decay curves of an ensemble of QDs in 3D photonic crys- 
tals with a continuous distribution of decay rates. We 
relate this distribution to the fact that QDs in various 
positions in the unit cell with random dipole orienta- 
tions experience different LDOSs. It is gratifying that 
recent calculations for Bragg onion resonators also con- 
sider nonexponential decay for similar reasonsi^l 

Our results demonstrate that large inhibitions and en- 
hancements of the spontaneous emission can be achieved 
with properly positioned and oriented efficient dipolar 
light sources inside 3D photonic crystals, at room tem- 
peratures and in large volumes limited only by the crystal 
size. The complementary case of a single QD in an opal 
has been recently studied in Refl28l. and for a single ori- 
ented quantum dot in a two-dimensional slab, interesting 
steps have been discussed in Refs. IT3lll5ll29l 
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